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o , 

. We discuss studies on determination of fragmentation functions and an application to 

^jQ' exotic-hadron search by using characteristic differences between favored and disfavored func- 

' tions. The optimum fragmentation functions are determined for pion, kaon, and proton in 

the leading order (LO) and next-to-leading order (NLO) of the running coupling constant 
Qs by global analyses of hadron-production data in electron-positron annihilation. Vari- 
' ous parametrization results are much different in disfavored-quark and gluon fragmentation 

' functions; however, we show that they are within uncertainties of the determined functions 

by using the Hessian method for uncertainty estimation. We find that the uncertainties are 
especially large in the disfavored-quark and gluon fragmentation functions. NLO improve- 
i-J^ ' ments are explicitly shown in the determination by comparing uncertainties of the LO and 

NLO functions. Next, we propose to use differences between favored and disfavored fragmen- 
tation functions for determining internal quark configurations of exotic hadrons. We make 
a global analysis for /o(980) for finding its internal configuration; however, uncertainties are 



(-H ' too large to specify the structure at this stage. 

§1. Introduction 

> 

Semi-inclusive hadron-production processes are important for investigating prop- 
00 ■ erties of quark-hadron matters in heavy-ion colhsions and for finding the origin 

of the nucleon spin in lepton-nucleon scattering and polarized proton-proton col- 



00 ' lisions. Fragmentation functions (FFs) are key quantities in describing such hadron- 



o 
o 



production processes in high-energy reactions. They indicate hadron-production 
probabilities from partons. 

The FFs have been determined mainly by hadron-production data of e^e~ re- 
, , - action. Many measurements were done in the Z'^-mass region, whereas lower-energy 

\ data are not sufficient. The determination of the FFs is not in an excellent situation 

d ' in comparison with the one of parton distribution functions (PDFs) , which is obvious 

from the fact that there are huge differences between the FFs of different analysis 
groups, especially for disfavored-quark and gluon FFs.M It led us to investigate un- 
certainties of the FFJ21i as it was done in the PDFs for the nucleon and nuclei.'^'''^ 
In addition, making global analyses in the leading order (LO) of the running cou- 
pling constant as and the next-to-leading order (NLO) at the same time, we showed 
that the role of NLO terms for reducing the uncertainties of the determined FFs. 
We provided a useful code in calculating the optimum FFs for the pion, kaon, and 
proton from our global analyses.'2l''Eli After our studies, there are works on related 
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global analyses of the FF^''^' and also on a hadron-model estimated 

Next, we proposed a possible method for exotic-hadron search by using the FFs!^ 
In particular, the FFs are usually classified into favored and disfavored functions. 
The favored means the fragmentation from a quark which exists in a hadron /i as a 
constituent in a naive quark model. The disfavored means the fragmentation from a 
sea quark. Therefore, internal quark configuration should be reflected in both FFs. 
This fact led us to investigate an interesting suggestion to use the FFs for a possible 
exotic-hadron search by looking at differences between the favored and disfavored 
functions.'^ 

We explain these works in this article. In Sec. [21 the FFs are defined in e^e~ 
annihilation processes. Our global analysis method is explained in Sec. El and results 
are shown in Sec. jH The idea of using the FFs for exotic-hadron search is introduced 
in Sec. [5l Our studies are summarized in Sec. [H 

§2. Fragmentation functions in e+e" annihilation 

The cross section of hadron-/i production in the e^e~ annihilation is described 
by a gg-pair production e+e" — )• qq followed by a fragmentation process from q {q or 
gluon emitted from q or q) to the hadron h. The FF is defined by the cross section 

F\z,Q') = ^'^-^^^^^^-^^, (2-1) 
(Jtot dz 

where utot is the total hadronic cross section, and is the virtual photon or 
momentum squared in e~^e~ — t- 7, Z^. It is equal to the center-of-mass energy squared 
s (= Q^). The variable z is defined by the energy fraction: 

- (2-2) 



^/i/2 Q ' 

where is the hadron energy. Namely, z is the hadron energy scaled to the beam 
energy {y/s/2). The fragmentation is described by the summation of all the parton 
contributions: 

F\z,Q^) = Y,Ci{z,as) D'l{z,Q^). (2-3) 

i 

Here, D^{z,Q'^) is a fragmentation function of the hadron h from a parton i (= 
g, u, d, s, ■ ■ •), Ci{z,as) is a coefficient function which is calculated in per- 
turbative QCD,'^''^^ and the convolution integral is defined by f{z) g{z) = 
dy/y f{y)g{z/y). 

The measurements of the FFs have been done in various Q^, whereas they are 
parametrized at a fixed point (= Qq) as explained in the next section. The initial 
functions at Qq are evolved to the experimental Q^ points by the standard DGLAP 
evolution equations. The equations are essentially the same as the ones for the PDFs 
by exchanging the splitting functions Pqg and Pgg. There are also some differences 
between their NLO expressions,-^ for example, in the modified minimal subtraction 
(MS) scheme. 
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§3. Global analysis method 

The FFs are expressed in terms of a number of parameters, whicii are determined 
by a analysis of the e"*" + e~ — t- /i + X data. The initial functions are provided at 
Qq as 

D^{z,Ql) = Ntz<'{l-zf', (3-1) 

where Nj^, a^, and (3^ are the parameters. An apparent constraint for the parameters 
is the energy sum rule: 

Y^Mt^Yl t dzzD^{z,Q^) = l, (3-2) 
h h 

where M/^ is the second moment of D^^{z,Q^). However, it is almost impossible to 
confirm this sum since the summation over all the hadrons cannot be taken prac- 
tically. In our analysis, we tried to be careful that the sum does not significantly 
exceed 1 even within analyzed hadrons 
There are two types in the FFs: fa- 
vored and disfavored functions. For the 
FFs of light quarks (u, d, s), they are as- 
sumed to be equal if they are favored or 
disfavored ones. The favored functions 
are given by 

Dl\z,Ql) = Df{z,Ql) 

= Nf z<^ {I - zf^-" , (3-3) 

for vr^. The vr^ productions from u, d, 
s, and s are disfavored processes so that 
they are assumed to be equal at Qq: 

Dt{z,Ql) = Df{z,Ql) 

= Df{z,Ql) = Df{z,Ql) 

= z"^^^ {I - z f^^ . (3-4) 

The FFs from a gluon and heavy quarks 
are defined separately as 

Dfiz,Ql) = Nfz''^^\l-zff, 

Df{z,ml) = Df{z,ml) = Nfz<\l - z)^^\ 

Df{z,ml) = Df{z,ml) = ivf - zff, (3-5) 

where rric and mf, are charm- and bottom-quark masses. The parameters in Eqs. 
(|3-3p . ()3-4p . and ()3-5p are determined so as to fit the data in Table HI where exper- 
imental collaborations, center-of-mass energies, and the numbers of data are listed 
for the charged-pion production. It is clear that most data are taken at the mass. 



Table I. Experiments, center-of-mass energies, 
and numbers of data points are listed for used 
data sets of e+ + — > 7r='= + X?* 



Experiment 


V~s (GcV) 


# of data 


TASSO 


12,14,22,30,34,44 


29 


TPC 


29 


18 


HRS 


29 


2 


TOPAZ 


58 


4 


SLD 


91.28 


29 


SLD (u,d,s) 


91.28 


29 


SLD (c) 


91.28 


29 


SLD (b) 


91.28 


29 


ALEPH 


91.2 


22 


OPAL 


91.2 


22 


DELPHI 


91.2 


17 


DELPHI (u,d,s) 


91.2 


17 


DELPHI (b) 


91.2 


17 


Total 




264 
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The FFs for the kaon and proton are parametrized in the similar way by con- 
sidering favored and disfavored functions. The amounts of data are almost the same 
as the ones in Table HI The detailed should be found in the original article.'^ The 
parameters for the kaon and proton are determined in separate analyses. 

One of our major purposes is to show the uncertainties of the FFs as explained in 
Sec. [TJ The uncertainties have been already estimated in the studies of nucleonic and 
nuclear PDFs.E^'® The same Hessian method is used for the uncertainty estimation. 
The is expanded around the minimum point ^: Z\x^(C) = X^(C + ~X^{0 = 
j Hij5(,i5(,j , where Hij is called Hessian which is the second derivative matrix, ^ 

indicates a parameter set, and ^ is the set at the minimum x^ point. The confidence 
region is given in the parameter space by supplying a value of Ax^- Using the 
Hessian matrix obtained in a x^ analysis, we estimated the uncertainty of the FF by 
[(5L>^^(z)]2 = Ax^ J2j,k [9D'liz)/dCj]^ HJ^ [dD'^{z)/dik\^ ■ There are some variations 
among groups on the appropriate Ax^ value for showing the uncertainty range in a 
global analysis. The details are explained in our articl^ about our Ax^ choice. 

§4. Determined fragmentation functions for pion, kaon, and proton 

We show obtained FFs of the pion by the 
X^ analyses of the e^e'^ — t- ti^X data in Fig. 
[U where the FF data in the form of Eq. (j2Tp 
and our parametrization result with an uncer- 
tainty band are shown at = M^. The good 
agreement with the data indicates that the fit 
is successful from small- to large-z regions. 

Next, each FF is shown for the pion on 
the left-hand-side of Fig. [5] with uncertainty 
bands in both LO and NLO (MS). It should 
be noted that the charm- and bottom-quark 
FFs are shown at the scale of their mass 
thresholds Q'^=m'^ or m^, whereas the others 
are shown at Q^=l GeV^. The uncertainties 
are generally larger in the LO, which indicates 
the FF determination is improved due to NLO 
terms. We notice that disfavored-quark and 
gluon FFs have large uncertainties in both LO and NLO. 

The NLO functions are then compared with other parametrizations of KKP, 
Kretzer, AKK, and DSS in Fig. [21 Some disfavored-quark and gluon functions, for 
example s-quark functions of Kretzer and AKK, are completely different between 
the analysis groups; however, they are within our uncertainty bands. There are 
not much differences between the groups in favored- (u), charm-, and bottom-quark 
functions except for the small-z region. The large discrepancies among the different 
parametrizations especially in the disfavored-quark and gluon functions were not 
clearly understood before our work. Therefore, it is important to point out in our 
studies that they are not due to an inappropriate analysis of some groups and that 
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Fig. 1. Charged-pion fragmentation func- 
tion in e'^e^ annihilation process. 
Our analysis results are compared 
with experimental data.'^ 
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Fig. 2. Determined fragmentation functions for the pion and their comparison with other 
parametrizations!^''^ The shaded bands indicate estimated uncertainties. 
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Fig. 3. Determined fragmentation functions for the kaon and proton 



they come from experimental errors and inaccurate flavor decomposition for light 
quarks. It is simply impossible to determine accurate disfavored-quark and gluon 
FFs from current experimental data. 

For the kaon and proton, our LO and NLO FFs are shown in Fig. [3l We 
found that both FFs are not better determined than the pion FFs. There is a 
similar tendency with the pion case that the disfavored-quark and gluon FFs are 
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not determined well. The NLO improvement, namely the uncertainty reduction, is 
clear in the kaon, whereas it is not apparent in the proton. We also found-^ in both 
kaon and proton that all the different analyses are consistent with each other in 
spite of their large differences in some functions because they are roughly within our 
uncertainty bands. 



§5. Exotic-hadron search by fragmentation functions 



Table II. Second moments of the FFs in the NLO 
for 7r+, K+, and p at Q^^l GeV^ 



Type 


Function 


2nd moment 


Favored 




0.401±0.052 


Disfavored 




0.094±0.029 


Favored 


or 


0.0740±0.0268 


Favored 




0.0878±0.0506 


Disfavored 


jjK+ 


0.0255±0.0173 


Favored 




0.0732±0.0113 


Disfavored 




0.0084±0.0057 



From the analyses of ordinary 
hadrons, pion, kaon, and proton, 
we found characteristic differences be- 
tween the favored and disfavored FFs. 
In Table [III the second moments are 
shown for vr"*", , and p. It is clear 
that the moments are larger for the 
favored functions than the ones for 
the disfavored functions. It indicates 
that internal quark configuration can 
be found by looking at flavor depen- 
dence of the FFs. This fact suggested 
us to use the FFs for exotic hadron searct 

The exotic means a hadron with internal quark configuration other than ordinary 
qq and qqq. This topic has been investigated for a long time; however, an undoubted 
evidence has not been found yet. In the last several years, there have been reports on 
exotic candidates mainly from Belle and BaBar collaborations in charmed hadrons. 
In our work, we investigated a possibility that the FFs can be used for an exotic 
hadron search by using differences between the favored and disfavored FFs. 

As one of the exotic mesons, we investigated a possibility of determining quark 
configuration of /o(980), which structure has been controversial for many years. Ac- 
cording to a simple quark model, it is described by the configuration [uu + dd) / \f2. 
However, it is known that theoretical strong decay widths are an order of magni- 
tude larger than the experimental widthjl-J Therefore, it is considered to be ss, 
tetraquark, or KK molecule state. It used to be considered as a glueball candi- 
date, but recent lattice QCD estimates indicate that the lowest scalar-meson mass 
is about 1700 MeV. There is an indication that the internal configuration could be 
determined by the radiative decay — ?■ /o7 and 27 decay /o — )• 27;^''^^' however, 
the FF method could become a better way for judging its internal structure as well 
as other exotic-meson configurations. 

We summarized in Table Hill how to judge the structure of the /o meson, espe- 
cially by the second moments and z-dependent functional forms. All the possible 
configurations are considered in the table. For example, if /o is an ss state, ss is 
formed from s by creating an ss pair from a radiated gluon. In the same way, a 
color neutral ss can be formed from a gluon by its splitting into an ss pair and a 
subsequent gluon radiation for color neutrality (see Ref.9 for details). The gluon can 
be radiated from either s or s, so that Mg could be larger than Mg as indicated in 
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Table III. Possible /o(980) configurations and their features in FFs at small 



Type 


Configuration 


Second moments 


Peak positions 


Nonstrange qq 


(uu + dd)/y/2 


Ms < Mu < Mg 


max ^ max 
^ Zg 


Strange qq 


ss 


Mu < Ms < Mg 


max ^ max 
Zu ^ Zs 


Tetraquark (or KK) 


{uuss + ddss) 1 \/2 


Mu^ Ms< Mg 


max max 
Zu ^ Zs 


Glueball 


99 


Mu Ms <Mg 


max max 
Zu Zs 



the table. If s is a favored quark, a significant portion of its energy (namely large z) 
is transferred to /o(ss). It leads to a functional form which is mainly distributed in 
the large-z region. In the table, this is denoted as z^""^ < zf"""^ . In the same way, 
the relations in the second moments and the functional forms are listed for other 
configurations. Our suggestions are intended to give an idea on the criteria, and 
further details need to be investigated in more sophisticated hadron models. 

There are data on the FFs of /o in the e"'"e~ annihilation. We have done a 
global analysis for determining the FFs in the same way with the analyses in Sec. 
[SP' This is intended to judge the quark configuration of /o by using the criteria 
of Table lllli The determined second moments are given by M„ = 0.0012 it 0.0107, 
Ms = 0.0027 ± 0.0183, and Mg = 0.0090 ± 0.0046. Then, the moment ratio becomes 
Mu/Mg = 0.43 lb 6.73. From the ratio 0.43, the /o seems to be mainly an ss state. 
However, it is obvious from its huge error 6.73 that a clear determination is currently 
not possible. 

The obtained FFs from the global analysis 
are shown in Fig. SI The up- and strange-quark 
functions are distributed mainly in the large-z 
region, which may indicate that u and s could 
be important constituents of /q. However, the 
uncertainties of the determined FFs are huge 
and they are an order of magnitude larger than 
the FFs themselves. Therefore, it is not possible 
to draw a conclusion on its structure from the 
global analysis at this stage. Hopefully, much 
better data will be reported in the near future 
possibly from the Belle collaboratioiPSJ for de- 
termining structure of exotic hadrons including 
/o(980). 
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Fig. 4. Determined FFs for /o(980).l^ 
Note that uncertainty bands are 
much larger than the FFs. 



§6. Summary 



Our studies on the fragmentation functions were reported. First, global analyses 
have been done for pion, kaon, and proton for determining their FFs. In the past, 
it was not clear why there are large discrepancies among various parametrizations 
on disfavored-quark and gluon FFs. We clarified that they are consistent with each 
other by estimating the uncertainties of the FFs by the Hessian method in the sense 
that all the distributions are within our error bands. We also clarified the role of 
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NLO terms in reducing the uncertainties by comparing the determined FFs and their 
uncertainties in the LO and NLO. Our code for calculating the FFs was supplied on 
our web site .'3' 

Next, we investigated a possibility of using the FFs for finding internal structure 
of exotic hadrons by using differences between the favored and disfavored FFs. We 
proposed to use the second moments and z-dependent functional forms of the FFs for 
determining the internal structure. As an example, the /o(980) meson was studied 
by considering all the possibilities of qq^ tetraquark, and glueball configurations. A 
global analysis has been done also for the FFs of /o by using the current e+e~ data; 
however, they were not accurately determined to the level of finding the internal 
structure, namely a difference between the favored-quark and disfavored-quark (and 
gluon) FFs. Much more accurate data are needed to discuss the internal structure 
by the FFs. 
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